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Abstract 

Qualitative information about breather initial profiles in the weak coupling limit 
of a chain of identical one-dimensional anharmonic oscillators is found by studying the 
linearized equations of motion at a one-site breather. In particular, information is found 
about how the breather initial profile depends on its period T. Numerical work shows 
two different kinds of breathers to exist and to occur in alternating T-bands. Genericity 
of certain aspects of the observed behaviour is proved. 

1 Introduction 

Chains of coupled anharmonic oscillators have many applications in condensed matter and 
biophysics as simple one-dimensional models of crystals or biomolecules. Of particular 
interest are the so-called "breather" solutions supported by such chains, that is, oscillatory 
solutions which are periodic in time and exponentially localized in space [Q. The simplest 
possible class of models, where all the oscillators are identical, with one degree of freedom, 
and nearest neighbours are coupled by identical Hooke's law springs has been widely studied 
in this context. The equation of motion for the position q n (t) of the n-th oscillator (n G Z) 
is 

Qn - a(q n+ i - 2q n + g n _i) + V(q n ) = (1) 

where a is the spring constant and V is the anharmonic substrate potential, which we 
choose to normalize so that V(0) = and V"(0) = 1. MacKay and Aubry have proved 
the existence of breathers in the weak coupling (small a) regime of this system [||. They 
noted that in the limit a — ► 0, system (|l]) supports a one-site breather, call it q = (q n )nez, 
where one site (re = say) oscillates with period T while all the others remain stationary at 
0. Using an implicit function theorem argument they proved existence of a constant period 
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continuation q a of this breather away from a = provided a remains sufficiently small and 
T £ 2vrZ+. 

Having established the existence of these breathers, the question remains: what do they 
look like? This question may be addressed by means of numerical analysis in several ways. 
One technique is to seek period T solutions of a spatially truncated (i.e. \n\ < N < oo) 
version of ([!]) by searching for fixed points of the period T Poincare return map using, for 
example, a Newton- Raphson algorithm ||. This method potentially allows the construc- 
tion of breathers far from a = 0, and can be used to determine the domain of existence 
of continued one-site breathers in the (a, T) parameter space (as applied by the authors 
to a non-standard discrete sine-Gordon system in H). This technique has its drawbacks, 
however. It is rather computationally expensive and there are subtleties concerning conver- 
gence of the Newton- Raphson method close to a = |j. On the other hand, if the weak 
coupling regime is of primary interest, then much useful information may be obtained from 



/ <9q Q 



(2) 

a=0 



the tangent vector at a = to the continuation curve Q a . This vector may be calculated 
by solving a certain initial value problem for a set of 4 coupled second order ODEs, a 
numerically trivial task. The small a behaviour of breathers can then be determined by 
approximating the curve q a by its tangent line at a = 0, that is 

q Q = qo + aqo + o(a). (3) 

Since this method is so computationally cheap, it is easy to make a systematic study of 
the period dependence of breather initial profiles for a variety of substrate potentials. This 
paper presents such a study. 

The paper is organized as follows. In section ^ we give a precise statement of the 
MacKay-Aubry existence theorem and reduce the calculation of q to numerically tractable 
form. In section || we present numerically generated graphs of the components of q against 
period T for various choices of potential, and extract from them information about breather 
initial profiles in the small a regime. The results show a generic behaviour of alternating 
T-bands of two qualitatively different types of breather, which we call "in-phase breathers" 
and "anti-phase breathers" . In section || we prove that this alternating behaviour is generic 
in a precise sense, and examine the small T limit analytically. Section || contains some 
concluding remarks. 



2 The direction of continuation of one-site breathers 

In the following we will assume that V : K — »■ R is twice continuously differentiable and has 
a normalized stable equilibrium point at (V'(0) = 0, V"(0) = 1). Consider the equation 
of motion for a particle moving in such a potential, 

x + V'(x) = (4) 

with x(0) = 0. Provided \x(0)\ is small enough, x(t) must be a periodic oscillation. All the 
potentials we consider will be anharmonic with classical spectrum (27r,oo). Anharmonic 
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means that the period of oscillation varies nondegenerately with x(0) > 0. The classical 
spectrum of V is the set of periods of the oscillations supported by V. 

Let xx(t) denote the solution of (Q) with period T > 2tt and x(0) > which has even 
time-reversal symmetry, xt{— t) = xj-(t). From this we may construct a 1-site breather 
solution of system (|l]) with a = 0, call it qo: 

1n,o(t) = { n7 , . (5) 

The MacKay-Aubry theorem establishes the existence of a continuation q a of this solution 
away from a = in a suitable function space, defined as follows. 

Definition 1 For any n € N, let C^ n denote the space of n times continuously differen- 
tiable mappings R — » K which are T -periodic and have even time reversal symmetry. Note 
that C^ n is a Banach space when equipped with the uniform C n norm: 

\q\ n := S up{\q(t)\,\q(t)\,...,\q ( - n \t)\}. 

tew 



Definition 2 For any n G N, 

:= {q : Z — > Cy n such that ||q|| n < oo} 

where 

||q|| n := sup |g m | n . 
iVote i/iai (^y n , || * ||n) *s a ^ so a Banach space. 

The required function space is J7^ 2 . By construction, qo € fiy 2 an< ^ * s exponentially 
spatially localized. 

Theorem 3 (MacKay-Aubry) If T £ 2-7rZ there exists e > suc/j i/iai for all a € [0, e) 
there is a unique continuous family q a 6 °^ solutions of system ([Jj at coupling a 

with qo as defined in These solutions are exponentially localized in space and the map 
[0, e) — > Q^ 2 c/uien i|/ a w q Q is C 1 . 

The idea of the proof is to define a C 1 mapping : fX£ 2 ® ^ ~~ * ^to> 

F(q, a) m = # m - a(q m+1 - 2q m + g m _i) + V'(q m ), (6) 

so that F(q,a) = if and only if q is an even T-periodic solution of system (|l|). In 
particular, -F(qo, 0) = by construction. Using T ^ 2ttZ + and anharmonicity of V, one can 
show that the partial derivative of F with respect to q at (qo,0), DF qo : £1^2 ~~ * ^ro» * s 
invertible (DF qo is injective with (DF qo )~ l bounded). Hence the implicit function theorem 
H applies and local existence and uniqueness of the C 1 family q Q satisfying 

F(q a ,a)=0 (7) 
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are assured. Persistence of exponential localization is proved as a separate step. 

The object of interest in this paper is qg € f2j 2 , the tangent vector to the curve q a 
at a = 0. This may be constructed by implicit differentiation of (|7j) with respect to a at 
a = 0: 



DF qo + 



dF 
da 







(q ,0) 
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dF 
da 



(qo,o) 



Straightforward calculation shows that 

[DF qo X ]m = 

for any x € Vt^ 2 , and 

(dF ^ 

(qo,0). 



Xm ~r Xm 
Xm + T/"(x T )Xr, 



m/0 
m = 




(8) 



0) 



(10) 



So evaluating the right hand side of (|8|) to find q^ = X G ^ 2 * s equivalent to solving the 
following infinite decoupled set of ODEs for {x m £ C^2 : m £ ^} : 



Xo + V"(x T )xo = -2x T 

Xm + Xm = XT \m\ = 1 

Xm + Xm = |m| > 1. 



(11) 

(12) 
(13) 

for \m\ > 1, so 



Recalling that T £ 2vrZ+ ( SO COS ^ C/jt 2 

) one sees from ( |i~3[ ) that Xm 
one need only compute xo and Xi (clearly x-i = Xi)- 

In each case we must find the even T-periodic solution of an inhomogeneous 2nd order 
linear ODE. For m = 0, 1 let y m (t) be the solution of the corresponding homogeneous 
equation with initial data y m (0) = 1, y m (0) = 0> so 



j/o + y"(x T )yo = o 
yi + vx = 0. 



(14) 
(15) 



(Of course y\ = cos.) Similarly, let z m (t) denote the particular integral of the inhomogeneous 
equation with initial data z m (0) = z m (0) = 0, so 



z + V"(x T )z = -2x T 

Z\ + Z\ = XT- 



(16) 
(17) 



Since Xmif) is even, it follows that Xm(t) = Xm(0)y m (t) + z m (t), and Xm(0) may be deter- 
mined by applying either of the T-periodicity constraints 



or 



Xm(T)= X m(0) Xm(0) 
Xm{T) = Xm(0) Xm(0) 



Z m (T) 

i - y m (T) 

z m {T) 

y m (T) ' 



(18) 
(19) 
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In either case, one sees that the tangent vector q , and in particular its initial value q (0) = 
(. . . , 0, 0, Xi(0),xq(0), xi(0), 0, 0, . . .) can be constructed once the 4 initial value problems 
(|T3-|17|) are solved on [0,T]. 



Having computed the initial value of the tangent vector, we may approximate the 
breather initial profile for small a using (^): 

(x T (0) + axo(0) + o(a) m = 
axi(0) + o(a) \m\ = l (20) 
o(a) \m\ > 1. 

So to first order in a, the continuation leaves all but the central (m = 0) and off-central 
(m = ±1) sites at the equilibrium position. The qualitative shape of the breather initial 
profile depends crucially on xi(0) (but not Xo(0))- If Xi(0) > 0, the continuation displaces 
the off-central sites from equilibrium in the same direction as the central site. The result is a 
hump shaped breather in which the central and off-central sites oscillate, roughly speaking, 
in phase (that is they attain their maxima and minima simultaneously). We shall call such 
breathers "in-phase breathers (IPBs)." If Xi(0) < 0, on the other hand, the off-central sites 
are displaced in the opposite direction from the central site, resulting in a sombrero shaped 
initial profile. In this case, the central and off-central sites oscillate, roughly speaking, in 
anti-phase (the central site attains its maximum when the off-central sites attain their 
minima, and vice- versa). We shall call such breathers "anti-phase breathers (APBs)." 
Sievers and Takeno ||] have argued (without giving a rigorous mathematical proof) that 
breathers of the latter type are supported by certain oscillator chains with no substrate 
potential and anharmonic nearest neighbour coupling, but their existence in systems of 
type (|I|) does not seem to have attracted attention previously in the literature. This is 
surprising in light of the numerical results described in section |3|. We shall find that for a 
generic substrate potential V, system (||) supports both IPBs and APBs, the type varying 
in bands as T increases through (2ir, oo). 



3 Numerical results 

Calculation of the constants Xo(0) and Xi(0) may be performed numerically by solving the 
initial value problems (|l4|), (|i"6| ) and ([D]) using some approximate ODE solver (recall the 
exact solution of (|l5|) is known). The results presented in this section were generated using 
a 4th order Runge-Kutta method with fixed time step 5t = 0.01. Each ODE requires the 
periodic oscillation xx{t) as input. For most potentials this function must itself be generated 
by numerical solution of the oscillator equation (||), which might just as well be solved in 
parallel with ( |i^|l6||T7| ) , making 4 coupled 2nd order ODEs in all. Since the correct initial 
displacement xt(0) for a given period T is not known in advance, the initial value problem 
is parametrized by initial displacement, T being determined from the approximate solution 
of @) by counting sign changes of xx{t) and linear interpolation. 

The results obtained depend crucially on whether the substrate potential has reflexion 
symmetry about the equilibrium position. 
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3.1 Asymmetric potentials 



The following asymmetric potentials were investigated: 



Morse: Vm(x) 
Lennard- Jones: Vlj{x) 
Cubic: V c (x) 



1 



2 
1 

72 



(1 



e~ x f 



12 



1 2 1 3 

2 3 



(21) 
(22) 
(23) 



Note that the equilibrium position for Vlj(x) is x = 1 rather than x = 0, as we have 
been assuming so far. Hence, even though Vlj is even, it is still asymmetric, since it is 
not reflexion symmetric about x = 1. Graphs of %o(0) and xi(0) against T for all these 
potentials are presented in figures 1 and 2 respectively. 




Lennard-Jones 



1 1 ' 1 1 1 

2 4 6 8 10 12 14 

T/ir 

Figure 1: Graphs o/xo(0) against T for various asymmetric substrate potentials. 

Figure 1 illustrates a clear difference in the large T behaviour of xo(0) between Vc and 
the other two potentials Vm and Vlj- Namely Xo(0) remains bounded as T — > oo for Vc, 
but grows unbounded for Vm and Vlj- This phenomenon appears to be linked to the large 
T behaviour of the T periodic oscillations xx{t). Due to the unstable equilibrium point of 
Vo at x = 1, xr(t) tends to a bounded homoclinic orbit as T — > oo for this potential, and 
Xo(0) is correspondingly bounded for large T. For Vm and Vlj no such bounded homoclinic 
orbit exists: xt(0) is an unbounded function of T and Xo(0) is correspondingly unbounded 
as T — > oo. In all cases, we note that Xo(0) > 0, so the continuation to a > displaces the 
central site further from equilibrium. 
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Figure 2: Graphs of Xi(0) against T for various asymmetric substrate potentials (solid: 
Morse; dashed: Lennard-Jones; dotted: cubic). 



From figure 2 we see that the sign of Xi(0) changes as T increases, leading to the 
prediction of T-bands of IPBs and APBs, as explained in section ^. The sign changes are 
clearly associated with the vertical asymptotes of the graphs at each T £ 27rZ + (note that 
these do not conflict with the contents of section |2| since the existence proof for q Q , and 
hence for q , breaks down when T € 27rZ + ). The presence of such asymptotes may be 
understood by using Green's function techniques to write down z\ (t) , the solution of ([H 



zi(t) = / sin(i — s)xt(s) ds. 
Jo 

Combining fl24| ) with (|l9|), and recalling that y\ = cos, yields a formula for xi(0)> 

fT 



Xi (0) 



i r 

— / cos(T — t)xT(t) dt. 

smTJo 



(24) 



(25) 



Equation @ 
Xi(0), unless 



shows that there is a vertical asymptote at T = 2nir, with a sign change in 



2nir 



COSt X2nn(t) dt = 0, 



(26) 



that is, unless X2nn(t) has vanishing n-th Fourier coefficient. It might appear from (^fj) 
that there should also be asymptotes at T = (2n + l)ir. Of course, this cannot be true 
since standard results on continuity of solutions of ODEs with respect to initial data imply 
that Xi(0) is continuous for T ^ 27rZ + . In fact, a simple argument using periodicity and 
evenness of xt demonstrates that 



(2n+l)7r 



costX( 2n+ i) 7r (t) dt = 



(27) 
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for all n G Z + and V. By contrast, we will prove in section |4| that (^) almost never holds 
(in a sense which will be made precise), so that the resonant periods T = 2nir generically 
separate IPB bands from APB bands. 

3.2 Symmetric potentials 

The following symmetric potentials were investigated: 

Frenkel-Kontorova: Vfk( x ) = 1 — cosrc (28) 
Quartic: Vq(x) = ^x 2 - ^x 4 (29) 
Gaussian: V G {x) = 1 - e~ x2/2 (30) 

The graphs in figure 3 show Xo(0) against T and closely resemble those of the previous 
section: xo(0) is bounded as T — > oo if xt tends to a bounded homoclinic orbit (Vfk an d 
Vq), and unbounded if xt(0) — ► oo as T — > oo (Vg)- 




Frenkel-Kontorova 



Quartic 



2 4 6 8 10 12 14 

T/lt 

Figure 3: Graphs of xo(0) against T for various symmetric substrate potentials. 

On the other hand, the graphs of xi(0) against T (see figure 4) are strikingly different 
from the asymmetric case. In each case, vertical asymptotes are present at T = 2nir only if 
n is an even positive integer. If n is odd, no asymptote is present. To explain this, note that 
evenness of V implies that xxit) is T/2 antiperiodic, that is Xx(t — T/2) = —xr(t), which 
in turn guarantees that equation (|26| ) holds whenever n is odd. The limit limr^nTr Xi(0) 
(n odd) exists and may, in principle, be computed using L'HospitaPs rule with 

— z x {T) = x T (T) + / [- sin(T - t)x T (t) + cos(T - t)f'(T)y (t)] dt, (31) 
where yo is the solution of (O) as before and f(T) := xr{0). 
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2 4 6 8 10 12 14 

T/u 



Figure 4 : Graphs of Xi(0) against T for various symmetric substrate potentials (solid: 
Frenkel-Kontorova; dashed: quartic; dotted: Gaussian). 



3.3 Almost symmetric potentials 

For m = 1 equation (|l8| ) becomes 



Xi(0) = T ^V (32) 
1 — cos 1 

From this equation it is clear that the behaviour of Xi(0) as a function of T is largely 
determined by the zeroes of z\(T). In particular, the sign of Z\(T) determines the sign 
of Xi(0). Figure 5 shows typical zero distributions for both asymmetric and symmetric 
potentials. It is interesting to note that no infinitesimal perturbation of the symmetric type 
distribution yields the asymmetric type distribution. Hence the dependence of Xi(0) on 
T for an almost symmetric potential must differ qualitatively from the typical asymmetric 
behaviour observed in section [3.1| . 

Asymmetric 

— ^ >p< 50[< >p< >o^ x-^ x-x- x-^ — - TV 7i 

2 4 6 8 10 12 14 16 18 



Symmetric 

— ^ ^ x-^< m x— ^ ^ x-^< ^ x— x- — *- T/7i 

2 4 6 8 10 12 14 16 18 

X = simple zero, 53 = double zero 

Figure 5: Typical zero distributions of zi(T) for asymmetric and symmetric substrates. 
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This observation motivates us to consider z±(T) for the following 1-parameter family of 
potentials: 

V e [x) = \x 2 - ix 4 - e -x\ (33) 

The zeroes of z\(T) in a region of the (T, log e) plane are plotted in figure 6. In the symmetric 
limit (loge — > — oo), the typical symmetric type distribution is recovered, as expected. For 
a fixed but very small e > (e.g. loge = —30), the small T oscillations xt(^) are pointwise 
essentially unchanged from the e = case, and the distribution for small T is numerically 
indistinguishable from the symmetric case. Only for large T does the oscillation xy(t) detect 
the asymmetry of V e , and a transition from symmetric to asymmetric behaviour occurs in 
the distribution. As e is increased, the period at which this transition occurs becomes 
smaller and smaller, so that the transition "propagates" leftwards as seen in figure 6. For 
sufficiently large e (e.g. loge > —3) the transition is complete and the typical asymmetric 
distribution is recovered. These observations appear to be essentially independent of the 
asymmetric perturbation considered. 




4 Genericity of numerical results 

The aim of this section is to explain some of the generic features observed in the numerical 
results of section |3| In particular we will consider the T —> 2tt + behaviour of Xi(0), and the 
generic presence of asymptotes at T S 2irZ + . 



4.1 The small period limit 

The numerical data suggest that Xi(0) — > +oo as T — > 2ir + for every potential V. We may 
confirm this prediction, at least in the case where V is analytic (as is the case in all our 
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examples), by means of a Linstedt expansion for xx{t) J?]]. We seek to construct a solution of 
the oscillator equation (|j) with initial data x(0) = e > 0, x(0) = 0, by power series expansion 
in e. Such a solution will be periodic in t with unknown period T(e) = 2ir/uj(e), where T(e) 
and hence uj(e) depends analytically on e. Defining a rescaled time variable s := uj{e)t 
(so the solution is s-periodic with period 2ir) and dependent variable y(s) := x{t)/e (so 
y(0) = 1, y'(0) = for all e), we substitute the expansions 

u(e) = a + aie + a 2 e 2 + . . . (34) 
y(s) = y (s) + eyi(s) + e 2 y 2 (s) + . . . (35) 

into (||) and Taylor expand V(x) about x = 0. Grouping terms by order in e, this results in 
an infinite set of coupled ODEs for {y n : n £ N}, to be solved in turn (n = 0, 1, 2, . . .) with 
initial data yo(0) = 1, yo(0) = y n (0) = Vn(0) = 0. Demanding that each y n (s) have period 
2tt then fixes the constants ao,ai,a2> • • •• For our purposes, it is sufficient to proceed only 
up to order e 2 , whereupon one finds that 

w(e) = l-Ae 2 + 0{e 3 ) 

T(e) = 2^[l + ^e 2 + 0(e 3 )], (36) 

where A = —V^(0)/16 — (V^ (0)) 2 /12 is positive for each of our examples (since they are 
all softly anharmonic, i.e. T > 2ir). Then 

x T{t) {t) = ecosw(e)i + e 2 yi (uj{e)t) + e 3 y 2 (^(e)t) + . . . (37) 

where all we need know about yi, y 2 , ■ ■ ■ is that they are even with period 2tt. We can now 
extract the asymptotic form of zi(T(e)) for small e: 

zi(T(e)) = i sm(T(e)-t)x T{e) (t)dt 
Jo 

1 f 27r . /2vr 



f 2n /2ir-s\ . 2 . , 

/ sin — -— ) (ecos s + e yi(s) + . . .) ds 
Jo V w(e) y 



Hence from equation (|32] 



W (£) 

A7r 2 e 3 + 0(e 4 ). (38) 



Xl(0) = T ^#T = ^ + O( e °). (39) 
1 — cos i (e 27Tj4e 



As e — > + , equations (pq) and (|39| ) reproduce the T — > 27r + asymptotic behaviour observed 
in all the numerical experiments. 



4.2 Vertical asymptotes 

We have seen that a sign change in xi(0); an d hence a transition from IPBs to APBs, must 
occur at T = 2mr (n € Z + ) unless zi(2mr) = 0, equation (|2"6|). In this section we will prove 
that, for a given n, the set of potentials on which (|2^) holds is negligibly small, so that the 
presence of asymptotes in the graphs of Xi(0) against T really is generic. To be precise, 
we will show that, for potentials in a certain Banach space V, the subset of potentials on 
which ( ]26| ) holds is locally a codimension 1 submanifold. 
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Definition 4 Let V = {V : R -> R suc/i i/taf V is C 2 and | \V\ \ 2 < 00} and F : C% 2 P -> 
6e i/ie mapping 

F(q,V)=q + V'(q). 
Note that F is a C 1 mapping between Banach spaces. 

Definition 5 A potential V £ V is anharmonic at q G C^ 2 if Q is nonconstant, F(q,V) = 
0, and ker DFr q y\ = {0}, where DFi q y\ : C^ 2 ~^ Cto denotes the partial derivative of F 
at{q,V). 

Clearly F(q, V) = means that q is a solution of Newton's equation for motion in poten- 
tial V. That injectivity of DF^ q y^ is equivalent to the standard definition of anharmonicity 
in Hamiltonian mechanics is shown, for example, in the original proof of Theorem |3] For 
all the potentials we considered, V is anharmonic at q = xt G C*t 2 ^ or an ^ ^ 7r - 

Definition 6 A perturbation neighbourhood of (q, V) G 2 ®F is a pair (U, f) where 
U C V is an open set containing V, and f : Li — > C^ 2 * s a C 1 map satisfying (a) every 
W G U is anharmonic at f(W) and (b) f(V) = q. Note that hi is trivially a Banach 
manifold. 

Lemma 7 For any V G V anharmonic atq G 2 i/iere exists a perturbation neighbourhood 
(U, f) of (q, V). The function f is unique on sufficiently small U. 

Proof: Anharmonicity implies that DFr q y\ is injective, and hence the standard solvability 
criterion in linear ODE theory guarantees it is also surjective (see section 3.3.2 of j|] for 
details). The open mapping theorem then ensures boundedness of DF,~L. Hence DF^ q> v) 
is invertible, and the implicit function theorem applied to F ensures existence of (U, f) and 
local uniqueness of /. □ 

Definition 8 For T G 2itZ + , we shall say that V G V is degenerate at q G C^ 2 if V is 

anharmonic at q and 

T 

cos t q(t) dt = 0. 

Theorem 9 Let T G 2nZ + , Vq G V be degenerate at qo G C^ 2 , an d (U, f) be a perturbation 
neighbourhood of(qo,Vo). The subset oflA on which degeneracy persists is a codimension 1 
submanifold. 

Proof: Consider the C 1 mapping / : U — ► R defined by 

I(V) = f T cost [f(V)] (t)dt 
Jo 

(L is C 1 since / is C 1 ). By the Regular Value Theorem (see appendix) the result follows 
if we establish that is a regular value of /, or in other words that I is a submersion at 
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every V € I 1 (0). Since ker Dly is of finite codimension in TylA = V it splits and hence it 
suffices to show that Dly : V — > R is surjective for all V G / _1 (0). 

Let V G I _1 (0)> 9 = /(V). For any SV G P, let Sg^y = -[DFT^JSV o g)] G C+ 2 . In 
other words Jg^y is the unique even, T-periodic C 2 solution of 

SQSV + V"{q(t))8 qsv = -SV\q(t)), 

whose existence follows from anharmonicity of V at q (invertibility of DF^ q y^ on Cy ). 
With this notation, 

DIy(SV) = / cost5q$y{t)dt 
Jo 

and surjectivity of Dly will follow if we exhibit SV G V such that DIy(SV) 7^ 0. We 
construct such a SV in two steps. 

For some to £ (0,T/2) and e > 0, let 6 G C^ 2 be a non-negative function with supp6 D 
[O.T/2] = [t - e, t + e]. Clearly 

f T ,T/2 

/ 6(t) cos tdt = 2 / 6(t) cos tdt^O 
Jo Jo 

(where T G 2irZ, + has been used) provided we choose to ^ cos _1 (0) and e sufficiently small. 
Since the critical points of q(t) are isolated, we may also assume that q is invertible on 
[to — e, to + e] with C 2 inverse q~ l : [x\ , 22] — > [*o — e, *o + e ] • By its construction the function 
g defined by 

J -[[DF {qy) b} oq-^x) xe[x 1 ,x 2 ] 
9{) \ x£[xi,x 2 ] 

is in C°(R) (since DF maps to C^ ). This allows us to construct a function 5V(x) G C 1 (R) 
such that Sq sv = b as follows. For x G ran(q), SV(x) = f Q x g(z)dz defines a C 1 function on 
ran(g). Since 5q sv is clearly independent of the behaviour of 5V outside ran(g), we extend 
SV to a C 1 function on R with support contained in some compact set K D ran(g). 

We cannot immediately conclude that Dly is surjective, since SV is C 1 but not nec- 
essarily C 2 . However, using a density argument it is straightforward to find a SV G C 2 
close enough to <5F such that Dly(SV) ^ stills holds. More precisely, Dly has a natural 
continuous extension, call it Dly to C . By construction SV satisfies DIy(SV) ^ and 
belongs to C\ (the C 1 functions with support contained in K). Since C\ is dense in C^-, 
there exists a sequence SV m G C\ C converging (in C 1 norm) to and by continuity 
DIy(SV m ) for all m sufficiently large. □ 

In the above, we have chosen V = (C 2 (R),|| • || 2 ) as our space of potentials. This 
choice was made for the sake of clarity and notational simplicity - many other choices 
would work. In fact, all but two (Vfk and Vg) of the example potentials considered in 
section || lie, strictly speaking, outside V, since they are unbounded. However, the analysis 
above can easily be adapted to deal with this: one simply replaces V by the affine space 
Vy := {W : \\W-V\\ 2 < 00}. 
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5 Concluding remarks 



In this paper, breather initial profiles in the weak coupling regime of a simple class of 
oscillator networks have been examined, focusing on the dependence on breather period 
T. The direction of continuation of one-site breathers was determined numerically using 
a simple and inexpensive numerical scheme. Two types of breather were identified, called 
IPBs and APBs. The numerical data suggest that generically these two types occur in 
alternating bands in the T parameter space (T £ (27r,oo)), and that the resonant periods 
T E 27rZ + separate an IPB band from an APB band. The genericity of this behaviour was 
proved rigorously. 

The distinction between IPBs and APBs, which does not appear to have attracted 
much attention in the literature, may have phenomenological implications in applications 
of the model (Q). One reason for interest in discrete breathers (particularly continued one- 
site breathers) is that, because of their strong spatial localization, they typically require 
little energy to achieve large amplitude oscillations close to the centre. This makes them 
good candidates for "seeds" of mechanical breakdown of the network, the idea being that 
the central oscillation becomes so violent as to break the chain (a similar mechanism is 
postulated as a mechanism for DNA denaturation, for example [||]). Since APBs stress 
the central intersite springs more than IPBs they presumably make more effective seeds of 
mechanical breakdown. 
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Appendix 

We recall some basic definitions and an elementary result of infinite dimensional differential 
topology, the Regular Value Theorem, which gives conditions under which the level set of 
a smooth function is guaranteed to be a manifold. 

Definition 10 We say that a closed subspace S of a complete topological vector space B 
splits if there exists another closed subspace C which is complementary to S, i.e. C + S = B 
and CnS = (0). 

Note that if B is a Hilbert space then any closed subspace S splits since we can take the 
complement to be C = S 1 - . Also if B is finite dimensional then since any subspace is closed 
all subspaces split. More generally any finite dimensional subspace (necessarily closed) of a 
Banach space splits, as does any closed subspace of finite codimension. 

Definition 11 A C 1 map f : X —>Y between Banach manifolds is a submersion at x G X 
if Df x : T X X — > Tj( x jY is surjective and ker Df x splits. A value y E Y is a regular value 
of f if f is a submersion for every x € f^ 1 (y). 
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Note that in the case that Y is finite dimensional then ker Df x is of finite codimension in 
T X X and hence is guaranteed to split. In this case we need only verify that Df x is surjective 
for / to be a submersion at x. 

Theorem 12 (The Regular Value Theorem) Let f : X — > Y be a C 1 map between Banach 
manifolds and y G Y be a regular value of f. Then f~ 1 (y) is a submanifold of X with 
T x f- 1 (y)^ker Df x . 

The proof follows from working in charts around x and f(x) and applying the Implicit 
Function Theorem. For details of a proof see M . 
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